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Abstract 

We discuss perturbative evolution of the polarized structure function g\ in the 
(x, Q 2 ) plane, with special regard to the small- x region. We determine g\ in terms 
of polarized quark and gluon distributions using coefficient functions to order a s . At 
small x gi then displays substantial scale dependence, which necessarily implies a corre- 
sponding scale dependence in the large x region. This scale dependence has significant 
consequences for the extraction of the first moment from the experimental data, reduc- 
ing its value while increasing the error. Conversely, the scale dependence may be used 
to constrain the size of the polarized gluon distribution. 
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1. Introduction 

Experimental data on the polarized proton structure function g\ [|l] have substantially 
improved within the last few months [HH, and results of comparable accuracy for the 
deuteron should be available soon [Q. Not only are statistical errors now substantially 
reduced 0, but also data which extend to rather small values of x (of order 0.003 0) are 
available, as well as determinations of the structure function at different scales for a given 
value of x. This suggests that it may be possible to see the effects of perturbative QCD 
evolution in the data. Indeed, the contribution of the polarized gluon distribution to g\ 
which appears in next-to- leading order, has been put forward || as one possible way of 
understanding the observed smallness of the first moment 

r!(Q 2 )= / 9l (x,Q 2 )dx; (1.1) 

JO 

it may now be possible to test for such a contribution by examining the x and Q 2 de- 
pendence of the structure function data. Also, the small- a; data @ indicate (albeit with 
large experimental uncertainty) a rise of the structure function at small x. Such a rise, 
even if not present in g\ at low scales, is generated at higher scales through perturbative 
QCD evolution || by a similar mechanism to that which drives the rise in the unpolarized 
structure function F^ at small x [[7] recently observed at HERA 

A quantitative understanding of the scale dependence and the small- x behaviour of 
the structure function is required in order to extract accurately its moments from the data: 
because data are taken at different values of Q 2 for each x bin, and over a limited range in 
x, it is necessary to evolve them to a common scale t = ln(Q 2 /A 2 ) and extrapolate them to 
all < x < 1 in order to determine any moment of g\ . The evolution is usually performed 
by the experimental collaborations 0|| by assuming that the virtual photon scattering 
asymmetry^ 

Al q 2 ) . = „(,, q») 2x P _ + y >1 (1.2) 

^1/2+^3/2 F 2 {X 1 Q 2 ) 

is scale independent, so that the scale dependence of g\ is given by that of the unpolarized 
structure function F 2 , which is known rather accurately. This can however be a rather 
poor approximation: as we shall see below, at small x d 1 ^ 1 diverges, and thus a large 



1 We follow the notations and conventions of ref. to which the reader is referred for theo- 
retical background. 
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error can be made if it is assumed that A\ is scale independent in that region. Likewise, 
the small- a; extrapolation is usually made by assuming Regge behaviour, which is not 
necessarily preserved by perturbative QCD evolution. This is especially dangerous since 
the contribution of the unmeasured small- a; region to the first moment Y\ eq. Ql.l|) is 
extrapolated from the lowest x data, which are necessarily taken at rather low values of 
Q 2 . The small- x contribution accounts for about 5% of the first moment |§, and because of 
large cancellations can give some 20 or 30% of the physically interesting singlet component 

ofiMl. 

A detailed investigation of the expected x and Q 2 dependence of g±, especially at 
small x, is thus called for. A first study of the scale dependence of the asymmetry eq. ( |L~2| ) 
was performed in ref. |I(J, where the one- loop Altarelli-Parisi evolution equations were 
linearized in hit to determine the approximate correction to be applied to the asymmetry 
data. The corresponding correction to the measured first moment turns out to be very 
small (of the order of a few percent, i.e. smaller than statistical uncertainties). However, 
the effect on individual data points may be appreciable; also, since these corrections are 
correlated, relatively small effects on single data points can lead to a significant change of 
the integral computed over a finite region of x. Of course, these effects would approximately 
cancel out if the integral could be measured over the whole range < x < 1 (due to the 
one-loop conservation of the first moment of g±). Such scale dependence corrections are 
thus increasingly significant as both the statistical quality and coverage of the small- x 
region of the data improves. It is then important to investigate the effects of the next- 
to-leading order coupling of polarized gluons to g±, in particular at small x, where the 
polarized gluon distribution will diverge. 

In the present paper we will analyze the x and Q 2 dependence of gi, including all the 
available theoretical [llj and experimental |||3| information. In particular, we will study 
the effects of including the coupling (which only appears at two loops) of the polarized gluon 
distribution to gi, and discuss the perturbative QCD prediction for the behaviour of g\ at 
small x. Even though a fully consistent two-loop determination of g\ is not yet possible due 
to the lack of knowledge of the full set of two-loop Altarelli-Parisi splitting functions, we 
will show that an approximate treatment suggests that evolution effects are significantly 
enhanced by the gluon contribution to g±, and that the inclusion of this contribution is 
necessary to understand current data. The size of these effects is essentially controlled 
by the size of the polarized gluon distribution; because this is unknown, it introduces 
an uncertainty in the determination of g\. However, especially in the small x region, 
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interference between the perturbative evolution of gluon and quark contributions may lead 
to large fluctuations of gi, so that the data may be used to constrain the size of this 
gluon contribution. We will discuss how present and future data can be used to determine 
the nature of the small- x behaviour of polarized parton distributions, thereby extracting 
information on the polarized gluon distribution. Also, we will discuss the impact of these 
evolution effects on the extractions from the data of the first moment of Ti, and show it 
to be not negligible at the level of accuracy of current experiments. 



2. Scale Dependence 



Our starting point is the relation between the polarized structure function gi and the 
polarized quark and gluon distributions Ag^ and Ag, which in general reads 1TT|JT3 



9i(x,t) 



(e 2 ) f'dy 



I T K S (f ' a s(t))AE(y : t) + Cf (J, a s (t))Aq™(y, t) 

+ 2n f C 9 (^a s (t))Ag(y,t) 



(2.1) 



where AE and Ag NS are respectively the singlet and nonsinglet quark distributions 



AE(x,t) = ^2{A qi (x,t)+Aqi{x,t)), Aq m (x,t) = £ 



' 4 - (e 2 ) 



i=i 



1 (e 2 ) 



(Aq i (x,t)+Aq i (x,t)), 
(2.2) 

and rif is the number of active flavours, with electric charge e^, (e 2 ) = Yl e 'i/ n f- The 
perturbative part of the x and t dependence of the polarized quark and gluon distributions 
is given by the Altarelli-Parisi equations [I3| : the singlet quark and the gluon mix according 
to 

a a (t) ^dy 



dt K ' ' 2% 



y 



P^,a s (t))AE(y, t) + 2n f P qg (^, a a (t))Ag(y, t) 



d . . a s (t) f dy 
—Ag(x,t) = — — / — 
dt 2tt J x y 



P OT (£ a s (t))AE(y, t) + P gg (* a 8 (t))Ag(y, t) 



(2.3) 



while the nonsinglet quark evolves independently as 



^Aq NS (x,t) 
dt H y 1 2tv 



a s (t) f x dy 



P^ q s (^a s (t))Aq^(y : t). 



(2-4) 



In eqs. (|2.1|) -( |272| ) we have taken both the factorization scale and the renormalization scale 
equal to Q , so that all scale dependence appears through t = In jfe. A different choice of 
factorization scale is discussed in ref. [11] . 
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The coefficient functions which relate the polarized parton distributions to g\ accord- 
ing to eq. (|2.1|) are given in the leading logarithmic approximation (used in ref. [10|j) by 



C^' s (x,a s ) = C^' NS (x,a s ) =5(l-x), C g (x,a s ) = 0, (2.5) 

so that g± is just a linear combination of polarized quark distributions, whose Q 2 depen- 
dence is entirely specified by eqs. (|2 . 3|) - ( PT^D (with the leading order splitting functions 
calculated in ref. |fL3f). 

Beyond leading order, splitting functions and coefficient functions are no longer uni- 
versal, hence even though the scale dependence of the (observable) structure function gi 
is determined uniquely, at least up to higher order corrections, its separation into con- 
tributions due to quarks and gluons is scheme dependent and thus essentially arbitrary. 
This ambiguity is somewhat constrained in the polarized case because of the anomalous 
conservation law satisfied by the singlet axial current ||14|| , whose nucleon matrix element 
is proportional to the singlet component of the first moment Fx eq. ( |1 . 1| ) of g\ ffl. The 
anomaly has the important implication || that the gluon contribution to g±, although 
formally of order a s , does not decouple as Q 2 — > oo. We thus expect the leading order 
approximation to g\ to be particularly poor in this case. 

Unfortunately, while all the coefficient functions have been calculated up to order 
a 2 in ref. [Oj, only the 0(a s ) corrections to the splitting functions P qg and P qq are 



known at present. Therefore we cannot yet perform a fully consistent next-to-leading order 
analysis. We will thus solve the Altarelli-Parisi eqs. ( |2.3| ) ( |2.4| ) at one loop, but include 
coefficient functions through order a s . The gluon coefficient function would be specified 
uniquely only once the next-to-leading contributions to the splitting functions are known 
in a particular scheme. We can partially fix it by insisting that the first moment of the 
polarized quark distribution be conserved by perturbative evolution. The first moment of 
the gluon coefficient function is then fixed by the anomaly equation to be JlJ 

J C g {x,a s )dx = (2.6) 

With this choice, the gluon density coincides with that which can be measured in inde- 
pendent hard processes |0| , in that all soft contributions to g± are absorbed in the parton 
distributions and do not contribute to coefficient functions [|IBJ. Furthermore, the first 
moment of AE coincides then with the nucleon matrix element of the canonical quark 



helicity operator [17]. We will further choose the higher moments of C„ to be equal to 



the first |fL8| . We may then test the sensitivity of our results to the inclusion of higher 
order corrections by choosing alternative forms of the gluon coefficient function which still 
satisfy eq. (|2.6| ), but result from different regularizations of infrared singularities |5,16 



In this respect our approach is the same as that adopted in ref. fL^| . As we shall see, 
evidence from the data supports the expectation that including the gluon coupling to g\ 
by enforcing eq. (|2.6| ) improves the quality of the leading order calculation, even though of 
course only a consistent two loop treatment will allow a definite conclusion. The inclusion 
of 0(a s ) corrections to the quark coefficient functions is known [pOlJlH to have a significant 
effect on the determination of Ti eq. ( |1 . 1| ) . 

With this choice, the Mellin transforms of the coefficient functions which we will use 
to determine g±, defined according to C(N, a s ) = f Q dx x N ~ 1 C(x, a s ) are 



Cf 5 (N, a s ) = l+^C F { {rf,(N) + 7) [3 + 2 (^(N + 2) + 7 )] 



An 



+2 



~6 



C 9 s (iV, a 8 ) = C NS (iV, a s ) + 0(a 2 s ) 
C g {N, a a ) 



47T 



2 ) 



(2.7) 



where ip(N) = dlogT(N) /dN . The explicit relation between the matrix element of the 
current and the first moments of the quark and gluon distributions is then 



(p,s\j»\ P: s) = Ms»a (Q 2 ); 



dx 



a. 



AX(x,Q A )-n f -^Ag(x,Q< 

Z7T 



where p, M and s are the momentum, mass and spin of the nucleon. The (scale dependent) 
coefficient of proportionality between ao and the singlet component of Ti can be determined 
in terms of the first moments of the quark and gluon coefficient functions comparing 
eq. 



with eq. (|2.1 



The anomalous dimensions analogously obtained by Mellin transform of the Altarelli- 



Parisi splitting functions in eq. (|2.3| ) are 

1 



ll q = C F 



3 

- + 



2 N(N + 1] 

99 _ r \— - V± 4- 4 

77V A L 6 9 iV(iV + 1) 



2(^(iV + l)+ 7 ) 



+ 0(a s 



2{ijj{N+ 1) + 7 ) 



+ 0(« s ), 



(2.9) 



where, for SU(3) color, Cp = I, Ca = 3 and Tp 



2 ' 



3. Small-x Behaviour 



We can now discuss the QCD predictions for the small- x behaviour of Aq, Ag, and 
thus gi. The treatment closely follows that of the unpolarized case ||. At small x, the 
Altarelli-Parisi evolution is dominated by the lowest moments of the splitting function |7| ; 
since these are singular at small enough N the leading behaviour is found by expanding 
the anomalous dimensions eq. ( |2.9|) in powers of iV around the rightmost singularity [^TJ : 



QQ 

W 


= C F 




+ 0(a s ) 


QQ 

Tn 


= T F 


-i + 2 + 0(JV) 


+ 0(a s ) 


l 9 N q 


= C F 




+ 0(a s ) 


1 9 n 9 


= C A 


"4 fn f 13 N 
. N V 9 6 / 


) +0(iV) + 0(e* s 



(3.1) 



Following the procedure described in ref. |§, i.e. diagonalizing the matrix of anomalous 
dimensions, inverting the Mellin transform, substituting the resulting splitting functions 



into the evolution equation (|2.3|) , and differentiating with respect to x, the eigenvectors 



(Aq^, A(7 ± ) are found to evolve in the (x, t) plane according to the pair of wave-like 



equations 



The relevant variables are 



In- 
to ' 



(3.2) 
(3.3) 



and xq and to = m (Qo/A 2 ) are arbitrary reference values of x and t, chosen in such a way 
that the approximate form of the splitting functions eq . ( |3TT|) is applicable for x xq and 
t ^ to- The coefficients y± and 5± are determined by the diagonalization of the singlet 
anomalous dimension matrix ( |3.1| ): the eigenvalues are X± = 2^o{j±/N + 5 ± + O(iV)), 
with 




(3.4) 



2(33-2n/), 

and Po = 11 — jTi/ is the leading coefficient of the QCD beta function. The eigenvectors 
v ^ = (Ag 1 * 1 , Ag^) are given (to leading order in N) by 



Ag ± = -2(lT^/l-^jA^ ± (3.5) 

Notice that the eigenvectors are not orthogonal. 

The small- x evolution equations satisfied by v ^ are thus of exactly the same form as 
those satisfied by G(x, Q 2 ) = xg(x, Q 2 ), where g is the unpolarized gluon distribution ||. 
There are two basic differences between the polarized and the unpolarized case. Firstly, 
since the polarized anomalous dimensions have their rightmost singularity at N = 0, while 
the unpolarized ones have it at N = 1, parton distributions here play the same role as 
parton distributions multiplied by x in the unpolarized case. Furthermore, since in the 
polarized case each entry in the matrix of anomalous dimensions has a simple pole at 
iV = 0, the polarized quark and gluon distributions evolve together on the same footing 
and the eigenvector of evolution is a linear combination of them, while in the unpolarized 
case, since only r y 9q and r y" are singular, the gluon distribution evolves independently 
according to a wave equation of the form ( |3.2| ), and the quark is then determined by it. 

These differences apart, the unpolarized and polarized cases have much in common 
since the form of the equation( |3.2|) is the same. In particular: 

(i) The equations ( |3.2[ ) are linear and causal. Their solution at (x ; , t') is determined by 
the knowledge of boundary conditions given on the x = xq and t = to axes for all 
x' < x < xq and t' > t > to. 

(ii) The equations are symmetric in £ and £, modulo the damping term proportional to 5±. 
Their solutions will display this symmetry, unless £ in which case the damping 
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term is increasingly important because the small- N expansion eq. (|2.9| ) is beginning 
to break down; any further asymmetry will be due to the boundary conditions, 
(iii) The general solution can be determined exactly in terms of the Green's function 
Io(2l±V£C) of the wave operator: 



(e,C) = Io(2T±VeC)e- 5± ^ ± (0,0) + J ^o(2 7 ±V / (^-e)C)e- 5±C ^ 7 ^ ± (r,0) 
+ | C rfC , /o(27±v / ^(C-C'))e 5±(C '- C) (^7^(0, O + ^(0, ('))■ 

(3.6) 

(iv) The asymptotic form of the solutions fall into two classes, according to the behaviour 
of the boundary conditions. Let us consider for definiteness the boundary conditions 
along the £ = axis, i.e., as a function of x at fixed t = to- For boundary conditions 
on both Aq and A g which are less singular than any power of x, the asymptotic 
behaviour for large £ and £ will be universal and given by the asymptotic form of the 
Bessel function Iq: 

(3.7) 

This behaviour holds as £C — > oo along any curve such that | — ► oo. If instead the 
boundary conditions have a power-like singularity of the form x~ x , the asymptotic 
result eq. ( |3.7[ ) only applies to the region £ < C v whereas for larger values of £ we 
have asymptotically 



v ± (^ C) ~ —=L== exp {2 7± v^C - S±C] ■ 



t; ± (£,C)~exp{A£+(4-5 ± )c}, (3.J 



i.e., the leading asymptotic behaviour reproduces the boundary condition. If the 
boundary conditions on Aq and Ag are different, but one (or both) has a power-like 
singularity, the most singular behaviour will always dominate asymptotically because 
of the mixing introduced by the evolution. 

The properties of the small-x solution to the evolution equations allow us to make some 
qualitative predictions on the expected behaviour of gi at small x. First, \gi\ is expected 
generically to rise at small x |J, at least as fast as eq. ( |3.7|) , i.e. stronger than any power 
of In - . This suggests that any rise of the boundary condition softer than this will soon be 
masked by the perturbative rise and thus not directly observable. Although both 7+ and 
7_ are positive, 7+ is larger so asymptotically the eigenvector v + eq. ( p75|) will dominate, 



and, at small x, Aq and Ag will eventually have the opposite sign. The sign of the rise of g\ 
is however not a priori determined, but rather is sensitive to the relative importance of the 
quark and gluon polarizations on the boundaries: if Ag is sufficiently large and positive, 
then asymptotically gi will tend to minus infinity. The onset of the asymptotic behaviour 
will lead, for some boundary conditions, to interesting interference effects. The shape of 
the rise will be determined by the qualitative form of the boundary conditions: it will take 
the universal form eq. (|3.7| ) for sufficiently soft boundary conditions, while it will tend 
to reproduce the rise in the boundary condition itself for harder ones. This is displayed 
in fig. 1, where we show a contour plot of the leading small- a; form of g\ computed by 
assuming that the Altarelli-Parisi equations reduce to their small-a; form eq. ( |3.2| ), with 
"soft" (fig. la) and "hard" (fig. lb) boundary conditions. If the boundary conditions 
are soft it is also possible to determine the asymptotic behaviour of the experimentally 
measured asymmetry eq. ( |1.2| ): given that the unpolarized structure function measured at 
HERA already displays the corresponding universal asymptotic form B], we have 



A x {x,Q 2 ) -xexp [(V2-l)2 7 VeC- yj^Cj , (3-9) 

where 7 and 8 are as given in ref. || , and we have neglected some very small n / dependent 
corrections. Thus, even though at small x the asymmetry falls as x, ln(-^4i) grows linearly 
with ^/ln(l/x) Int. 

Putting everything together, we expect that if sufficiently numerous and precise data 
on gi at small x and large Q 2 are available (say, x < 0.1 and Q 2 > 10 GeV 2 ) then it should 
be possible, just as in the unpolarized case ||, to determine experimentally whether the 
small- x structure function and parton distributions are soft or hard. Clearly, this is not 
yet the case, since only a couple of the SMC data points and none of the E143 data || 
fall within this range. This is demonstrated explicitly in fig. 1, where a scatter plot of the 
small x data is superimposed on the leading predicted asymptotic small x behaviour for 
soft and hard boundary conditions: it is clear that the data are insufficient to distinguish 
between the two scenarios simply because only a few of the experimental points fall close 
to the asymptotic small-a; region. 

We expect however that even with a few data at small x and reasonable Q 2 it should 
be possible to constrain the size of the gluon polarization, given the very strong sensitivity 
of g± and its scale dependence to it. Finally, because the small x eigenvector is a linear 
combination of quark and gluon (with a large coefficient multiplying the gluon contribu- 
tion), and because both the quark and the gluon distributions grow at small x but with 
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opposite signs, we expect in general the appearance in gi of fluctuations due to interfer- 
ence between quark and gluon contributions, until the asymptotic behaviour sets in. The 
observation of such fluctuations would put a strong constraint on the shape and size of the 
gluon, although presumably now in a non-universal way (i.e., dependent on the particular 
form of the boundary conditions). 



4. The Data 

In order to establish some more quantitative results, we need direct input from the 
experimental data. We use the SMC and E143 determinations of the scattering asymmetry 
eq. ( |1 . 2| ) , together with the NMC determination of the unpolarized structure function 
F2 J2!I and the SLAC determination |2^] of the ratio R(x, Q 2 ) of the longitudinal to 
transverse virtual photoabsorption cross section, to obtain values for g\(x,Q 2 ) according 
to eq.( |1.2| )i. This provides us with a determination of gi along a set of points (x^, Q 2 ) in 
the (x, Q 2 ) plane; due to the kinematic coverage of the present experiments these points lie 
for each experiment on a curve which is a monotonically rising function Qe xp (x). We then 
choose a physically motivated parametrization of the initial polarized quark and gluon 
distributions Aqi and Ag as a function of x for all < x < 1 at a reference starting 
scale Q 2 = Qq', gi is then determined for all Q 2 > Qq, and in particular at each data point 
(xfc, by means of eq. ( pTT| ) and the evolution equations (|2.3| ),( [2^|) . The free parameters 
in the initial parametrization are determined by a best fit of the evolved distribution to 
all the available data. Notice that, because the Altarelli-Parisi evolution is causal, the 
structure function determined in this way for all points (x,Q 2 ) such that Q 2 > Q 2 xp (x) 
is independent of the specific parametrization, to the extent that the data determine the 
shape of gi along the line Qexp( x )- The only real model dependence here comes in the way 
gi is decomposed into quark singlet and nonsinglet and gluon contributions. 

We parametrize the initial quark and gluon distributions as 

A£(x, Ql) = N (a q , q , a q ) r) q x a * (1 - xf" (1 + a q x) 
A(? NS (x, Ql) = N (a NS , /9ns, a NS ) ??nsx Qns (1 - x)^ NS (l + a NS x) (4.1) 
A#(x, Ql) = N (a g , p g , a g ) rjgx™ 9 (1 - xf 9 (1 + a g x), 



2 For consistency we do not include kinematic higher twist corrections in determining g\ from 
Ax, as was done in ref. Q, since these are not included in the extraction [f22| of F2 from the data, 
and because there exists no systematic treatment of higher twist corrections to either g\ or F2 ■ 
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where N(a, f3, a) is fixed by the normalization condition 

N(a,j3,a) f dx x a (l - xf{l + ax) = 1. (4.2) 
Jo 

Similar parametrizations have been used in recent global fits of polarized parton distri- 
butions ||19|| . The large x behaviour of the parametrization eq. Q4.1| ) is inspired by QCD 
counting rules from which we expect (3 q = /3ns — 3 and (3 g ~ 4. The polynomial 
in x (which we take to be first order) is a phenomenological interpolation which has the 
purpose of improving the quality of the fit in the intermediate x region. 

Several estimates for the small x behaviour of gi are available. Arguments based on 



the dominance of known Regge poles |25j lead to the expectation 



9i ~ x x ; < A < 0.5. (4.3) 

x— >0 



However, coherence arguments |24| suggest that (at a typical nucleon scale) the polarized 



gluon distribution should be related to the unpolarized one g(x) according to 

As 



~ x. (4.4) 

g x^o 



Thus, if the behaviour of is g dominated by a soft pomeron [26 j, then g(x) ~ - so that 



the lower bound for A in eq. ( fO| ) is saturated, but if g(x) has a harder behaviour then 
A < (A ~ —0.1 with a supercritical pomeron, and even lower with the Lipatov hard 
pomeron ||26|| , however values as low as —0.5 are now clearly ruled out by data from 
HERA |0]). A model of the pomeron based on nonperturbative gluon exchange gives 
the still singular but softer behaviour 

g x ~ -21nx. (4.5) 

x— >0 



Finally, it has been argued [^8[ that it is possible for negative signature cuts to induce an 



extremely singular behaviour 

9i ~. -rV- (4-6) 
xm x 

Nonperturbative arguments do not generally distinguish between the leading be- 
haviour of the quark and that of the gluon. We will hence consider several representative 



cases of soft or hard boundary conditions of the form eq. (|4.1| ) for both the quark and 
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gluon distributions.! Notice, however, that since QCD evolution mixes the quark and 
gluon distributions (as in eq. ( |3.5|) ), if either of the two has a very singular behaviour, on 
propagation the stronger singularity will dominate both. Also, behaviours of the form of 
eq. Q4.5|) or eq. (|4.3|) with A & 0.3 will grow rapidly upon evolution according to eq. ( |3T7| ) 
and thus their precise form will be hidden by the perturbatively generated rise of g±. 

We choose the starting scale at which the parametrization eq. ( |4.1| ) is imposed to be 
Qq = 1 GeV 2 . Here, this has the advantage that Qq < Ql xp (x) for all x, so that we may 
compare to the data without having to evolve g\ backwards in Q 2 (which would in general 
be unstable, due to the causal nature of the Altarelli-Parisi equations). Of course the price 
to be paid for this choice is that at such a low scale there could be sizable higher twist 
effects. However, in the unpolarized case these effects seem mostly to be concentrated in 
the large x region [p9| , where Q 2 xp (x) is large. Because the evolution of the leading twist 
is decoupled from that of the higher twist contributions, the overall effect on our results 
of having neglected higher twist corrections should be small. However, our final results 
for g±(x, Q 2 ) should only be trusted in the region where the higher twist effects are small, 
and in particular do not necessarily provide an accurate description of the data at large x 
(x <^ 0.3) and small scales (Q 2 ~ Qq ). 

At Q 2 — Qq, the quark distributions should be dominated by the valence component, 
hence we will take (3 q = /3ns an d a q — <2ns- We also take for simplicity the gluon parameters 
to be the same as the singlet ones, i.e. (3 g = (3 q and a g = a q , on the grounds that QCD 
evolution will anyway level small differences in the initial values of the various coefficients. 
More detailed fits could be obtained by relaxing these assumptions or increasing the order 
of the polynomial in x; however, our aim is not to propose a new parton parametrization, 
but rather to study the scale dependence of g±, especially in the small x region, where 
these details should make very little difference. The accuracy of present day data is in any 
case not good enough to justify fitting all these parameters independently JOJ. 

Given that no clear-cut theoretical prediction for the small x behaviour of Aq and Ag 
at a fixed scale exists, while we do not expect the present data to allow fitting the exponents 
a q and a g , we will let them vary independently over the range of values —1 < a q , a g < 0, 
and study the dependence of the results on this choice. The range of positive values 



We will not consider logarithmic corrections to the power-like behaviour of the form (4.5) 



or (4J3) since the coefficient functions already display a logarithmic rise at small x, and thus it 
would be very hard to disentangle from the data such corrections to the boundary conditions. 
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eq. (|4.3| ) all give essentially the same behaviour of g± (because of the logarithmic rise in 
the coefficient function) and are thus all accounted for by the flat case a q = 0, a g = 0. 
The lower bound is necessary since the first moment of the parton distributions must be 
finite. The nonsinglet quark distribution, which does not mix with the gluon and with 
the sea, will be assumed to be valence-like, i.e. to have «ns > 0; specifically we will take 
a NS = 0.2, which is the best-fit value of ref. |l(J (where only the earlier data in the valence 
region were considered). 

We are thus finally left with seven free parameters, namely the two parameters f3 and 
a controlling the large x behaviour, the two small-x exponents a q and a g , and the three 
normalization coefficients r\ qi and r\ g . The normalization of the nonsinglet, however, 
is fixed by the fact that the first moment of Ag NS is Q 2 -independent to all orders because 
of the conservation of the nonsinglet axial current. It can thus be determined from the 
knowledge of the nonsinglet component of the matrix element of the axial current measured 
in hyperon (3 decays: 



[ Aq m (x,Q 2 )dx=lg A + {a 8 , (4.7) 
Jo 



where [28 



a a = 1.2573 ±0.0028 

(4.8) 

a 8 = 0.579 ±0.025. 

The normalization of the gluon and singlet quark contributions to gi cannot be fixed 
analogously, due to the anomaly which afflicts the singlet axial curent. At a low scale, 
such as Q 2 ~ Qq, r\ q might be expected to be close to the normalization of the polarized 
quark distribution (i.e., the total fraction of spin carried by quarks) as calculated in quark 
models. One may estimate this from the normalization of the octet current as eq. (JOT), 
which for valence quarks according to the Zweig rule will coincide with that of the singlet. 
The difference between r\ q and the observed small value of Y\ (which is much larger than 
the evolution effects) is then be made up by a large value of r\ g || . This interpretation of 
the data, however appealing, need not necessarily be true; we will therefore consider three 
different scenarios: i) a "maximal" gluon case, namely r\ q = as; ii) the opposite limit where 
rj g = and any polarized gluons are generated dynamically; iii) a more general fit where 
both r/ q and r/ g are left as free parameters. Note that rather than using the experimental 
determinations of the first moment of g± as a constraint on r\ q and r\ gi as in ref. ||19|| . 



we 

will instead determine F± using our own fits to the g± data, evolved to a common scale. 
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The results of the fit to the data are displayed in table 1 and figure 2 (in several typical 
cases) and can now be used to draw some quantitative conclusions on the scale dependence 
and small- a; behaviour of g±. First, it is clear that a statistically significant evolution effect 
is seen in the data. Indeed, throughout the intermediate x region 0.03 < x < 0.2, where 
both experiments have data, there would appear to be a systematic discrepancy in that 
all of the E143 data fU lie systematically about one standard deviation below the SMC 
data || . The difference is however seen to be entirely accounted for by the fact that in each 
x bin the E143 data are taken at a smaller value of Q 2 than the SMC data. This is displayed 
explicitly in fig. 2a-c: in particular, fig. 2c shows that even if r\ g = and a q = 0, i.e. if the 
small x behaviour of the initial conditions is very soft, the dynamically generated gluon 
polarization is sufficient to drive the required amount of perturbative evolution. Notice 
that (as shown explicitly by the calculation of ref. ||10|| ) not enough evolution would be 
found if we had set C g = in eq. (|2.1j). In fact, the total x 2 of the fit shown in fig. 2c 
(which is the least sensitive to this effect since it has minimal gluon) increases by more 
than three units if we set the coefficient functions to their one loop values eq. ( |2.5| ). This 
provides a posteriori evidence in favor of our approach, i.e. it supports the assumption that 
the inclusion of the next-to-leading order coupling of Ag to g± even when the Altarelli- 
Parisi equations are only solved at leading order improves the determination of g\\ this 
can be thus interpreted, albeit with some caution, as direct experimental evidence for the 
anomaly-induced contribution of Ag to g\. 

Even though the evolution of Y\ is fixed in a universal way by the anomalous dimension 
of the axial current j§ (and is slight, since it starts at two loops), the evolution of gi for 
any finite value or range of x is a leading order (potentially large) effect. However, it is not 
universal, but rather depends on the specific form of Aq and Ag. In practice, the parameter 
which controls the exact amount of evolution is essentially the value of r) g , i.e., the size 
the gluon contribution at the starting scale. Large values of rj g lead to significantly larger 
evolution. The value of r\ g cannot be determined from the data, but may be somewhat 
constrained. The limit on rj g may be read off table 1; its actual value depends on the values 
of a q and a g . In general, the data favor intermediate values of r) g , i.e. such that r\ g > 0, 
but still not large enough that the condition rj q = ag be satisfied at the starting scale. 

On the contrary, no significant constraint is put by the data on the values of the 
small x exponents a q , a g , as expected from fig. 1. Thus a good fit can be obtained with 
either "hard" and "soft" boundary conditions, such as, for example a q = a g = —0.5 
and a q = a g = (see fig. 2a, b). In fact, with all combinations of small- x exponents 
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it is possible to accurately fit the data (table la) by suitably adjusting the value of r\ g . 
However, for different classes of values of a q and a g the qualitative form of the perturbative 
evolution in the intermediate x range (0.03 < x < 0.1) may be significantly different. In 
particular, if a g < a q (gluon more singular than quark) gi becomes rapidly large and 
negative at small x. It follows that only small values of r\ g (fig. 2a) are allowed so that 
the decrease of g\ falls outside the observed range. If instead a q < a g (quark at least as 
singular as gluon), then the onset of the small- x asymptotic behaviour (where AG and 
Aq are anticorrelated, according to eq. (|3j|)) tends to generate oscillatory interference 
patterns in g± (fig. 2d). These oscillations are more pronounced if r/ g is large (hence the 
evolution is stronger), and for soft boundary conditions (i.e. a q ~ a g <i —0.3), when the 
asymptotic double logarithmic growth eq. (|3.7|) only sets in slowly. This also allows us to 
rule out a maximal gluon for soft boundary conditions, since then fluctuations would be 
too large (fig. 2e). Finally, if a q is very large, then it is possible to satisfy the condition 
r\ q = ag even though in fact r\ g is very small (compatible with zero): in such case g± grows 
very large in the unobserved small- x region. This scenario |30] is also compatible with the 



data (fig. 2f): an arbitrarily large contribution to T± concentrated at very small x does 
not influence the evolution at larger x and thus cannot be seen in g\ (x) measured over any 
finite region of x. 

Conversely, if r/ g is fixed (for example on the basis of theoretical prejudice), then the 
values of a q , a g can be fitted (or at least constrained). In particular, if r\ g is very large 
(maximal), then only values a q ,a g < —0.5, with a q < a g are allowed, the corresponding 
minimum of the x 2 m the (a q , a g ) plane being rather narrow. If rj g = 0, then a valence-like 
quark is favored (a q ~ +0.5) , albeit with a sizable uncertainty. For intermediate values of 
r\ g the x 2 is very smooth in the (a q , a g ) plane, and a q and a g are essentially unconstrained. 

In order to check that these results are not a byproduct of our treatment of next-to- 
leading corrections (which are included in the coefficient functions but not in the Altarelli- 
Parisi splitting functions), we have repeated our fits with the two alternative forms of 
gluon coefficient functions given in ref. |HJ. The results do not change in a statistically 
significant way. 
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5. The First Moment 

On the basis of the quantitative knowledge on the evolution of g\ acquired in the 
previous section, we may now discuss the values of its first moment I\ eq. ( |1 . 1|) , and in 
particular the contribution to it from the small x region 

rf(Q 2 )= [°' 01 dx gi (x,Q 2 ). (5.1) 
Jo 

Values of Ti and rf x at various scales, and the value of the singlet contribution ao(Q 2 ) 
eq. (|2.8| ) are given in table 2, for several of the choices of the fit parameters which give 
good agreement with the data. 

The contribution eq. ( |5.1| ) of the small x extrapolation to Ti, at the scale of the 
data, turns out to be approximately the same for all reasonable fits (apart of course for 
the extremely singular case of fig. 2f). The uncertainty in the knowledge of the small-x 
behaviour of the starting quark and gluon distributions eq. fl4.1| ) is thus seen to have very 
little effect on the value of T±. 

However, there is a hitherto unnoticed sizable uncertainty in the value of T\ at any 
given scale which is only indirectly related to the small x extrapolation and is due to the 
scale dependence of g±. This may be seen by comparing, for instance, the results of fitting 
the data with rj g = and a q = (fig. 2c), or with a q = —0.5, a g = and maximal gluon, 
Vq = a 8 (fig- 2d). These fits have comparable x 2 , thus providing an equally accurate 
description of the data, and the same value of rf x (<5 2 ) at Q 2 = 3 GeV 2 . Nevertheless 
they lead to values of ri(Q 2 ) which differ by almost 10%, hence to values of its singlet 
component which differ by almost a factor 2. 

The origin of this uncertainty can be understood by remembering that whereas the 
singlet contribution ciq eq. ( |2.8| ) to Ti is scale independent at leading order, the integral of 
gi restricted to a limited x region is not. Inspection of the results of table 2 and fig. 2 shows 
that, for rjg > 0, in the intermediate x range (0.03 ~ x ~ 0.2) gi, and thus the contribution 
to Ti from this region, is a rather rapidly increasing function of Q 2 ; this is also manifestly 
displayed by the data. This implies that the contribution Tf^ from the (unmeasured) small 
x region must decrease as Q 2 increases. However, this Q 2 dependence of rf x is extremely 
sensitive to the value of r] g : if rj g is large, then rf x decreases very rapidly, due to the 
large oscillations induced in gi, while this is not the case if r\ g is small. But the value of 
rj g is only weakly constrained by the data, since the corresponding evolution effects are 
mostly concentrated in the experimentally unaccessible small x region. This implies that 
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if rj g is large, then Ti can be substantially smaller at the starting scale because g± in the 
intermediate x region, which provides the bulk of it, has been evolved more. 

As a corollary to this result, we find that our determinations of Fi are systematically 
smaller than those quoted in refs. [U|§. This is because the data points which provide 
the bulk of Ti are taken at a scale Q 2 larger than the nominal average scale of either 
experiment: for example in ref. || Q 2 xp (x) is always larger than 20 GeV 2 for x above 
0.1, while (Q 2 ) = 10 GeV 2 for this experiment. Since in this region g\ increases as the 
scale increases the approximate treatment of the scale dependence in refs. ||,[| leads to 
an overestimate of Ti.il Conversely, the small-x extrapolation of ref. @ overestimates the 
data, because it is obtained by extrapolating the smallest x data points, which are taken 
at very low Q 2 (~ 1 GeV 2 ), and then assumed to apply at < Q 2 >= 10 GeV 2 ; this turns 
out not to be the case for ref. || partly because their value of < Q 2 >= 3 GeV 2 is closer 
to that of the small x bins, and partly just because a straight line drawn from the smallest 
x points of this experiment happens to be closer to our results. 

This analysis is summarized pictorially in fig. 3, where we display the scale dependence 
of both Ti and the small x contribution to it, rf x , for representative minimal, intermediate, 
and maximal values of r\ g . Clearly, Tf x decreases significantly as a function of scale if r\ g 
is maximal, whereas it is essentially scale-independent if it is minimal (the slight rise 
displayed in the figure in such case is not significant to the level of accuracy of the present 
treatment: it could be turned into a decrease for different forms of the next-to-leading 
corrections which we do not include). As a consequence, larger values of r/ g correspond 
to smaller values of IV Since however the value of r\ g cannot be fixed on the basis of 
theoretical arguments, and is not constrained by the data, this is to be considered as an 
intrinsic uncertainty on the determination of V\. However, even with minimal values of 
inclusion of the next-to-leading evolution effects discussed here leads to a smaller value of 
Y\ than that found through the approximate treatment of the scale dependence given by 
the experimental collaborations 

4 Notice however that part of the effect is due to the fact that our determination of g± from 
the data of ref. || is systematically smaller than that presented there, because of the kinematic 
higher twist effects included in their determination of g± from the measured asymmetry. 
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6. Discussion 

The main result of our analysis is that the next-to-leading order coupling of the polar- 
ized gluon distribution to g\ has a sizable effect on the scale dependence of the structure 
function. This effect has important consequences for the extraction of the first moment 
Ti of <7i from the data because, even though T\ is scale-independent at leading order, the 
contribution to it from any finite x range (such as the experimentally accessible range) is 
not. In particular, the contribution to Ti from the unobserved small-x region turns out 
to be a decreasing function of the scale, this decrease being compensated by an approx- 
imately equal increase of the contribution from the measured region. The importance of 
the evolution effects is largely independent of the detailed shape of the quark and gluon 
distributions, but increases with the size of the polarized gluon distribution. 

As a consequence, due to the neglect of these evolution effects, the value of Tx is 
somewhat overestimated by present experiments, and the uncertainty, which is due to 
the lack of knowledge of the size of the polarized gluon, is substantially underestimated, 
more accurate approximate < Q 2 > closer to the value of Ql xp (x) in the x region gives 
the dominant contribution to Ti. Indeed, this considerably in reducing the uncertainty 
on the region. The largest uncertainty would then contribution which would be evolved 
uncertainty in the knowledge of the small x behaviour of polarized parton distributions 
leads to a comparatively much smaller ambiguity, at least excluding the possibility of a very 
large contribution to g± from the very small x region (such as a delta-like contribution). 
A more accurate determination of Ti and especially of the associated error can only be 
obtained by including the full scale dependence in the analysis of the data. From table 2 
we suggest the value 

T^Q 2 = 3 GeV 2 ) = 0.120 ± 0.005 ± 0.005 (6.1) 

where the first error is statistical and the second is an optimistic estimate of the theoret- 
ical uncertainty in the scale dependence due to the uncertainty in the size of the gluon 
distribution. The systematic error, which is typically of order 0.010, must still be added. 

The intrinsic ambiguity related to evolution could be substantially reduced if more 
precise data in the intermediate x region (0.03 ^ x ^ 0.2) were available for somewhat 
different values of Q 2 : because of the instability of gi under perturbative evolution even 
a moderate improvement in experimental knowledge could result in much more stringent 
theoretical constraints on the size of the polarized gluon distribution, which in turn would 
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be enough to determine the small x behaviour of g\ to much better accuracy. Such infor- 
mation would thus be more conclusive than direct information on gi in the small x region 
itself. The determination even of the qualitative behaviour of g\ at small x, in the same 
way as in the unpolarized case || would instead only be feasible with a kinematic coverage 
of the small x region comparable to that which HERA data provide at present for F2. 

In short, the scale dependence of g± induced by the axial anomaly turns out to be an 
essential ingredient in its phenomenology: a more detailed understanding of it will only 
be possible once the actual size of the gluon polarization is known more accurately. On 
the theoretical side, the determination of the full set of two loop Altarelli-Parisi splitting 
function will allow a fully consistent treatment. It will then perhaps be possible to extract 
information on the polarized gluon distribution, which is of considerable intrinsic theoret- 
ical interest, directly from structure function data even over a moderate range of x and 

Q 2 - 
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Table 1: Best-fit parameters and x 2 for the fit of g\ to the data of refs. [Qj3[: the letters 
a)-f) refer to the plots in fig. 2. 



20 







Ti(3) 


Ti(10) 


rf(i) 


rf (3) 


rf (io) 


a (3) 


Fitted gluon: 


a) 
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Table 2: Values of the first moment T 1 (Q 2 ) eq. ((T^), at Q 2 = 1, 3, 10 GeV 2 , its small 
x component eq. ( |5.1| ), and the matrix element ao of the singlet axial current for several 
typical fits from table 1: the letters a)-f) refer to the plots in fig. 2. 
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Figure Captions 



Fig. 1 Contour plot of In g\ computed from the small x evolution equation eq. in the 
plane of the variables defined in eq. (|3.3| ), with Qq = 1 GeV and xq = 0.1. The 
polarized gluon and singlet quark distributions are evolved up from the form eq. ( fO] ) 
with (3 = 3, r] q = 0.6, r\ g = 2.5 (maximal gluon), a = and a) a q = 0, a g = —0.2 
(soft); b) a q = -0.5, a g = -0.5 (hard). A scatter plot of the SMC (+) and E143 (x) 
data is superposed for reference. 

Fig. 2 Evolution of the structure function gi compared to the SMC (crosses) and E143 (di- 
amonds) data. The full set of parameters for the various fits are given in table 1. 
a) fitted rj g , hard boundary conditions; b) fitted rj g , soft boundary conditions; c) dy- 
namically generated gluon {j] g = 0); d) maximal gluon, hard boundary conditions; 
e) maximal gluon, soft boundary conditions; f) Ellis-Jaffe quark with small gluon 
contribution. 

Fig. 3 Scale dependence of Ti eq. (|1.1|) (a) and the the small x (x < 0.01) contribution to 
it eq. ( |5.1|) (b). Solid line: fitted gluon (as in fig. 2a); dot-dashed line: dynamically 
generated gluon (as in fig. 2c); dashed line: maximal gluon (as in fig. 2d). 
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